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ABSTRACT 

The hyperfine splittings Ac = {rriE,* — rnnj — {rnD*+ — mD+) and A^ = 
(msj — rriej — {niB'O — uibo) are analyzed in the framework of an effective 
lagrangian possessing chiral, heavy flavour and spin symmetries, explicitly bro- 
ken by a complete set of flrst order terms. Among these terms, those responsi- 
ble for the difference between the couplings gp*p*n and gp*pn are evaluated in 
the QCD sum rules approach. Their contribution to Ap, and to A^ appears to 
quantitatively balance previously estimated chiral effects nice agreement with 
the experimental data, solving a suspected puzzle for heavy quark theory. 



1 Introduction 

The spectroscopy of heavy mesons is among the simplest framework where the ideas and 
the methods of heavy quark expansion can be quantitatively tested. Recently, attention 
has been focused on the combinations fl], ||, ||, Q : 

Ad = {rriD* - m^J - (mz).+ - mD+) (1.1) 

Afi = (rriBi - niBj - (m^.o - niBo) (1.2) 

which are measured to be |^: 

Az3 ^ 1.0±1.8 Mel/ (1.3) 

Ab^1.0±2.7 MeV (1.4) 

The above hyperfine splitting is free from electromagnetic corrections and it vanishes 
separately in the SU{3) chiral limit and in the heavy quark limit. In the combined chiral 
and heavy quark expansion, the leading contribution is of order irts/mQ and one would 
expect the relation |jl|: 

Tfl 

Ab = —Ad (1.5) 

nib 

In the so called heavy meson effective theory [0, which combines the heavy quark ex- 
pansion and the chiral symmetry, there is only one lowest order operator contributing to 
Ad.b: 

A2O2 = ^Tri^X^i/^Km^i)!^ (1.6) 

where i, j are heavy flavour indices and a, 6 light flavour indices. The 4x4 Dirac matrix 
if* describes the spin doublet P, P*, with P heavy meson composed by the heavy quark 
Qi and the light antiquark g^. The matrix m^ is 

"^e = (^"^9^ + ^ W^) (1-7) 

Here mq is the light quarks mass matrix and ^ = exp{iM / f) , where M is the pseudoscalar 
3x3 matrix and / the pseudoscalar decay constant (we take / = 132 MeV) . 

By taking nis/AcsB — 0.15 and by taking A2 — ^qcd — O.lGeV^ one would estimate: 

,(2) 
D 

Af ~ 6 MeV (1.9) 

Given the present experimental accuracy, the above estimate is barely acceptable, as an 
order of magnitude, for A^, while it clearly fails to reproduce the data for A^. If the 
contribution from O2 were the only one responsible for the hyperfine splittings, agreement 
with the data clearly would require a rather small value for A2. 

In chiral perturbation theory, an independent contribution arises from one-loop correc- 
tions to the heavy meson self energies [Q, evaluated from an initial lagrangian containing. 



A^^ ^ 20 MeV {l.i 



at the lowest order, both chiral breaking and spin breaking terms. The loop corrections 
in turn depend on an arbitrary renormalization point /i^ (e.g. the t'Hooft mass of di- 
mensional regularization). This dependence is cancelled by the /i^ dependence of the 
counterterm A2(/i^)02, as it should happen for any physical result. A commonly accepted 
point of view is that the overall effect of adding the counterterm consists in replacing /i^ 
in the loop corrections with the physical scale relevant to the problem at hand, A^^g. 
Possible finite terms in the counterterm are supposed to be small compared to the large 
chiral logarithms. With this philosophy in mind, two classes of such corrections has 
been estimated in ref. and they give (for the values of the parameters given by these 
authors) : 

A^ ~ +30 Mel/, A]y ~ +65 MeV, (1.10) 

A^ ~ +10 MeV, A^ ^ +22 MeV, (1.11) 

Here A° represents the contribution of the chiral logarithm and A^ is a non analytic 
contribution, of order m^^^. This provides a rather uncomfortable situation since, to 
account for the observed data, one should require an accurate and innatural cancellation 
between (A'' + A^) and the finite terms from A^'^\ contrary to the usual expectation. 

The chiral computation giving A*^ + A^ is however incomplete 0, because it does not 
include the spin breaking effect due to the difference between the P*P*7r and the P*P7r 
couplings (P = D,B), defined by the relations: 

< vr^(g) P''{q2)\P*-{qi, e) > = 2 gp,p^ ^e^ • q, (1.12) 

In 

< n{q) P*(g2,e2)|P*(gi,ei) > = -i^gp*p*n e^.^^e^'e^g^gf (1.13) 

Jit 

The scaling law 2gp*pT^mp/ f^^ for the strong D*D7T coupling constant was first proposed 
in 1^, |]. The splitting between the coplings ( |1.12| ) and ( |1.13| ) is of order l/mq, and 



therefore has to be taken into account in the chiral computation, to work consistently at 
the desired order. 

In the present paper we will provide an estimate of gp*p*Tr — gp*p-K based on a QCD 
sum rule, and, by including this additional spin breaking effect, we will complete the 
evaluation of A^^^ coming from the chiral loops. 



2 The Hyperfine Splitting 

To better clarify the importance of gp*p*n — Qp'P-k for the problem at hand, we remind 
that the effective lagrangian for heavy mesons and light pseudoscalars, at first order in 
m'^ and in the light quark masses rUq reads: 

C = Cq + C'' + C^ (2.1) 



Here Cq represents the chiral, heavy flavour and spin symmetric term: 

£ = -iTr[Hlv^d^Ht] + ^—Tr[d^T)d^Y.] 

o 

+ \gTr[H:Hl^,^,]{Ani (2.2) 



A, = ed,i-id^e (2.3) 



where S = ^^ and: 



From the last term in eq. ( ^.21 ), one obtains the P*Ptt and P*P*tx couplings defined 
in eqs. ( p..l2| ) and ( |1.13|) , in the limit nip — > cxd: 

gp*p-K = gp*p*-K = g (2.4) 

The leading chiral breaking corrections are given by: 

a = AoTr[mgS + S"^mg] 

+ KTr[HlHt]{m^)l (2.5) 

The second term in eq. ( |2.5D is responsible for the mass splitting between strange and 
non-strange heavy mesons: 

A, = 2Aim, (2.6) 

One has approximately A^ ~ 100 MeV", Ai ~ 0.33. 

The third term, listed for completeness, gives an equal contribution to each heavy 
meson mass, it does not affect the hyperfine splitting, and it does not play any role in our 
analysis. 

Finally the terms of order I/tuq, breaking either the heavy fiavour or the spin sym- 
metries, are given by: 



^^ = -^Tr[H:a,.HX'']im-Q 



i^i 



+ f ^Tr[if:/f^7.75](mQ^)i(^-)^ 

+ '{^-^Tr[H:^,^,H^]{mQyMnt (2.7) 

The first term in eq. ( p.7|) is responsible for the splitting A between the 1~ and 0~ heavy 

meson masses: 

2A , , 

A = 2.8 

3 



For the B, B* system A ~ 46 MeV , whereas for D, D* A is approximately 141 MeV, so 
that one has: 

A~0.10-0.11 Gel/^ (2.9) 

The second term in ( |2.7|) breaks only the heavy flavour symmetry, making the B*B^^'n 
and D*D^*^TT couplings different. The third term breaks also the spin symmetry and 
contributes differently to the P*P7r and to the P*P*tt couplings. This is precisely the 
effect relevant to the hyperfine splitting. To this order in l/mq one has: 

gp*p*. = 9(l + —] 9P'Pn = 9(l + ) (2.10) 

and 

a — b 
Ag = gp*p*^ - gp*Pn = g (2-11) 

The chiral and spin symmetry breaking parameters relevant to the hyperfine splitting 
are the light pseudoscalar masses m.^, mx and ?7i^, A^, A and A^. In terms of these 
quantities, one finds P, 0, ffl: 



A; 



rt^A r A2 A2 A2 



_IA 






ml ' ■' ■ m^ ■ ■■ ■ ml 



+ tI^itiP'''""*-'^' 



The dependence upon the heavy flavour P = D, B is contained in the parameters A and 

A,. 

The first term in eq. ( |2.12| ) is the so called chiral logarithm 0] . In the ideal situation 
with pseudoscalar masses much smaller than Acsb, it would represent the dominant 
contribution to Ap. For the case of D and B mesons the corresponding values have been 
listed in eqs. ( [I.10| ) and ( |1.11|) as A^j and A^, respectively. There a value g^ = 0.5 has 
been used. 



The second term in eq. (|2.12|) represents a non analytic contribution of order m 



3/2 



which, although formally suppressed with respect to the leading one, is numerically more 
important, because of the large coefficient 247r. It is given by A\) and A]^ in eqs. (|1.10|) 



and (1.11 



Finally, the last term in eq. ( 2.12 ) ^ is also of order m^/^. It can be numerically 



important as soon as Ag/g is of order 10% and, if equipped with the right sign, it can 
cause a substantial cancellations of the previous two contributions. 



3 QCD Sum Rules for gp*pT^ and gp*p*T^ 



The coupling gp*p-„ has aheady been calculated in ^ by means of QCD sum rules and 
here we proceed to a similar computation concerning the coupling gp*p*T^. We start from 
the correlator: 

A^^Mi.q) = 1 I dx< 7r(g)|T(V;(x)V;t(0)|0 > e'^^^^ = A{qlqlq^)e,,^pq^q( + ... (3.1) 



where V^ = wy^Q is the interpolating vector current for the P* meson. 

We compute the scalar function A in the soft pion limit g — > 0. This implies gi = g2 
forcing to use a single Borel transformation, and it is the origin of the so called parasitic 
terms 0. The correlator in ( p.l|) can be calculated by an Operator Product Expansion: 
we keep all the operators with dimension up to five, arising from the expansion of the 
current V^{x) at the third order in power of x, and the heavy quark propagator to the 
second order. The result is: 



A{qlql^) 



U 



+ 



1 



ql - ml {qf - mlY 



2 < uu > nib Sf^ml 



3/. 



9 



+ 



(q! 



m 



2^3 



lOm^/^mf ml < uu > m^ 



3/. 



+ 



m,Q < uu > m,l 



{ql 



m, 



2\4 



u 



(3.2) 



In eqs.( p.2|) < mm > is the quark condensate (< uu >= — (240MeV^)^), mg and mi are 
defined by the equations 

< ugscr ■ Gu >= rriQ <uu> (3.3) 

< 7r{q)\uD^-f^-f5d\0 >= -iUmlq^ (3.4) 



and their numerical values are: ttiq = 0.8 GeV^, 



m^ 



0.2 GeV^ PI, ITT 



Proceeding in a standard way, we now compute the hadronic side of the sum rule. We 
can write down for A{q\, q^., 0) the following dispersion relation: 



A{ql,qlO) = —Jdsds'-^ 



p{s,s') 



q!)is' 



ql) 



(3.5) 



It should be observed that we have not written down in 



subtraction terms because, 



as proven in |0, only a subtraction polynomial P^{qi,q2) could be present in 
it would vanish after the Borel transform. 



but 



We divide the integration region in three parts 0. The first region (I) is the square 
given by ml < s,s' < sq and it contains only the B* pole, whose contribution is 



AiiqlqlQ) 



—2gB*B*TvfB*'^B* 

fAqI-^B')iqi-mB) 



(3.6) 



where fs* is defined by 

< 0\V^iO)\B*{e,p) >= ejB*mB^ (3.7) 

The second (II) integration region is defined as follows: ml < s < Sq and s' > Sq or 
iTT'b < -s' < So and s > sq. Here we obtain a contribution coupling the vector current V^ 
to the pion and the B*. Introducing the form factor V as 

< n{q)\V^\B*{q,,e) > = V{ql)e,,^pe' q^^ (3.8) 

where q2 = Qi — q, ^g get 



Auiql ql 0) = fs^ms^ { -^^^ + /^''^ , ] (3.9) 



^^ "^ Uj = jE'lfl-B* I ^ . . 

qt-niB*" q^-uiB* 



In the previous formula one does not have to include the B* pole contribution to ^(q'^), 
being already taken into account in Aj. We assume that, taken away B*, a single higher 
resonance of mass m' contribute to V 



Vresiq') = ^ (3.10) 



where k is an unknown constant. 



The third region is defined by s, s' > Sq, and under the assumption of duality it should 
coincide with the asymptotic limit g^ = g| ^ — oo in ( |3.2| ). One gets: 



Aui{qlqlQ) = ^/^- (3.11) 

The hadronic side of the sum rule is the sum of the contributions from the three 
regions: 

'^^'' ^'' ^ = /.(gf-mlO(gi-ml) + {q!-ms*'){q!-m'^) + J^ ^^'^^^ 

We have put gi = g2 and k' = mB*fB*k. 

Equating now the hadronic and the QCD sides of the sum rule, respectively given by 
eq. (|3.12| ) and ( p.2| ), and taking the Borel transform with parameter M^ we find: 

'^gB-B-Tv'mB* Ib* I ;/ I ( X/l\j'2\/f i l\ 

— ^^ + k +exp{-S/M ){f^-k) = 

1 (2 <uu> rrih SfT,m{\ 



exp iVt/M' 



u- 



M2 \ 3/^ 9 



1 f bmlf-j^ml rriQ < uu > mi\ 1 m^ < uu > m\ 



M4 V 9 QU J M6 6U 

expQ/M^S{M^) (3.13) 
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In the previous formula we have put m ~ sq and we have introduced the parameters 



mt 



5 = Sq — niB* and fi = m^. — v/t^. 

Differentiating ( 3.13|) respect to the variable 1/M^ and combining the first and second 
derivatives in order to eliminate the unknown parameter k', we obtain the following sum 
rule: 



gB*B*TT 



U exp Q/M^ 



2mB.2/B.2 



6 



+{2Q + 6)di,M2S{M^) + dln,2S{M^)\ (3.14) 



To eliminate the parameter k' one could also combine the first derivative with the 
original sum rule ( |3.13| ): we have checked that the two procedures give the same numerical 
results. We have used the second derivative to make an easy comparison with the sum 
rule for qb^bit 0: 



gB*Bii 



Af^mlrriB* 



exp (fi'/M^ 



X 



X 



2mlfBfB*{SmB^' + ml) {6' - 5'/\b'b/M^ - /\b*b) 
n'{n' + S')S'iM^) + (2fi' + 5')di/M'^S'iM^) + dl/M^S\M^)] (3.15) 

Q — Ab*b and 



where Ab*b 


= ms*^ - 


ml 


,S' = 


So - ml = 6 + Ab*b, fi' = m| - 


-ml = 




S\M^) 


= 


u- 


<uu> 






3mhf^ 






+ 


1 

M2 _ 


2 <uu > rrib lO/^rTif m,^ <uu> 
3U ' 9 ' 3m,U 






+ 


1 


lOmlf^^ml ml <uu > m^ 
9 ' QU _ 


+ 




2M4 






+ 


ml< 


uu > ml 





6/.M6 



(3.16) 



Eq. ( |3.15| ) differs slightly from the one given in 0, since it keeps track of the mass 
difference A^.^. The sum rules ( |3.14|) and ( p.l5|) have to be analyzed in the duality region, 
i.e. the region in M^ where there exists a hierarchy among the different contributions of 
higher dimension operators (this fixes the lower bound for M ^ ) ; moreover we impose that 
the contribution of the parasitic term does not exceed that of the resonance term, which 
fixes the upper bound in M^. In this way we obtain for the B M^ in the range 20—40 GeV^ 
(for So = 33 - 36 GeV^ ), and for the D M^ = 4 - 7 GeV^ (for sq = 6 - 8 GeV^ ). Using 
m^h = 4.6 GeV and m,c = 1.34 GeV one gets: 



p2 
Jb* 9B*B*tt 

Id* 90*0*71 
and for the gp*pn coupling 

Ib Ib* QB'Btv 
fo fo* gO*On 



0.0094 ± 0.0018 Ge1/2 
0.017 ± 0.004 Ge\/2 ; 

0.0074 ± 0.0014 GeV^ 
0.0112 ± 0.0030 Ge\/2 



(3.17) 



(3.18) 



Once multiplied by Impj Jt^ the figures in ( |3.18| ) agree with those given in 



A recent calculation of the quantity reported in Eq. ( ^.181 ) has been given in Ref. 
T3| . When expressed in our units their results are as follows: fs fs* gB*Bn = 0.0079 ± 



0.0007 GeV^ and fn Id- gD*Dn = 0.018 ± 0.002 GeV^. The result for the B is only 
slightly larger than our outcome Eq. ( p.l8| ), whereas the result for the D is significantly 
(^ 60%) larger. The origin of the discrepancy is in the different range of values for the 
Borel parameter M^, that, in the case of Ref. |]13|, are generally smaller. A possible origin 
of this difference is the fact that, while in this paper we use QCD sum rules in the soft pion 
limit, in I IB] light cone sum rules are adopted, which results in an expansion in operators 
of increasing twist instead of increasing dimension. In particular we have included a 
dimension 5 contribution which is proportional to the ml <uu> condensate. This term 



has no counterpart in |jT3|; since it has to be kept small, its inclusion in [T^ might result 
in a more stringent constraint on the hierarchy among the different contributions of the 
Operator Product Expansion, and, therefore, in a more stringent lower limit on M^. 



We now expand the sum rules ( |3.14| ) and ( |3.15| ) in the parameter l/rriQ, keeping the 
leading term and the first order corrections. The leading term is the one surviving in the 
limit nib — ^ C)0, and has already been calculated in for gp*pn- 

To extract the l/mq corrections we introduce the following parameters, finite in the 
large mass limit: 

(3.19) 



E 



M' 



yo 



So -ml 



u = mp — m^b 



2m.b 2m,h 

and the l/mg corrections to the leptonic decay constants, gp^p-,, and gp*p*n 



L 



M 




V^ 



I 



M* 



F ( A'\ 



(3.20) 



The coefficients A and A' have been computed in |]14| [0, but only for the B mesons. 
The corrections found are large and suffering of large uncertainties: there are significant 
numerical differences between |]14[ and ||15|| . 



In the limit r/ig -^ oo the right hand sides of ( |3.14| ) and ( ^.15] ) coincide (notice that 
S" = S" in this limit), confirming the result anticipated in (|2.4|) and giving: 



gF' 



f^ exp (uj/E) 



(yo 



UJ 



where 



So{E) = U 



yoSo{E) + (yo + co)^^/ESo{E)) + 9i%So(^) 
<uu> bfT,m.\ 



+ 



TTIq < uu > 



(3.21) 



(3.22) 



SUE 36^2 48^3/^ 

The duality region extends for E = A ~ 6 GeV and yo = lA — 1.3 GeV. Numerically one 
obtains: 



F^ g = 0.040 ± O.OObGeV^ 



(3.23) 



While this result agrees within the error with that given in |^, the central value 
reported here is 15 % larger due to a slightly different choices of the phenomenological 
parameters. 



We can then write the l/niQ expansion for the function S and S': 
SiE) = SoiE) + ^^ S'iE) = SoiE) + ^'^^^ 



where 



S,{E) 
S[{E) 



mq 






mq 



(3.24) 



9^ 



< uu > bf^mi uIq < uu > 



3U 9E ASE^U 

From (|3.14|) and ( p.l5|) , one gets the following sum rules for the parameters a and b: 

U exp uj/E 



(3.25) 



2gE\yo-uj) 



{iu' + AX){yoSo{E) + d,/ESoiE)) + 



+ 2yoUjS,{E) + 2(1/0 + uj)dyESi{E) + 2dl^j,Si{E) 
u? + 4A u;2 + 4A 



+ 



+ 



2E 



+ 



2(2/0 - ^] 



-2(A' + tu) 



(3.26) 



/^expcu/E 
2gF\y^-u) 



uj\yoSoiE) + d,/ESoiE))- 



+ 2yocoS[{E) + 2(2/0 + co)di/ES[iE) + 29i%S;(E) 
00^ + 4A 00^ + 4A 



+ 



+ 



+ 



{A + A' + Auj) 



(3.27) 



2E 2(2/0 - to) 

where A has been given in ( p.9|) . From the previous sum rules one gets: 

a + 2A' = -0.15±0.20GeV b + A' + A = ~1A5 ±0.20 GeV (3.28) 

and 

a-b+{A' -A) = 0.99±0.02GeV (3.29) 

Notice that the difference has a quite smaller uncertainty, due to a partial cancellation of 
terms depending on the threshold. 

Neglecting radiative corrections, A and A' are given by P 



A = -uj + ^ + 3G^ 



A' = -- + ^-G^ 
3 2 



(3.30) 



Notice that the splitting of the couplings depends on the quantity a — b that contains 
only the difference A' — A given by: 



A' -A = -iJ- 4Gs 
3 



(3.31) 



There is disagreement in the hterature on the values of the parameter Gy,'. at the h 



quark mass scale from ref.[|T4[ one gets Gs = (0.042 ± 0.034 ± 0.023 ± 0.030) GeV , while 
in ref. |T^ the central value Gs ~ —(0.052) GeV is quoted. In view of this discrepancy, 
to provide an estimate of the difference ( p.31| ), we will approximate A' — A k. 2/3a; k, 
0.4 GeV , obtaining 

a- 6^ 0.6 GeV (3.32) 



4 Discussion and conclusions 



From ( p.ll| ), (|3.32|) and from the formula ( p.l2|) of the hyperfine mass splitting we obtain: 



Ab ^ r(27.3 + 61.4 - 75.8) MeV = 12. V MeV (4.1) 



Notice that we have used in eq. ( p.l2|) f = f^ = 132 MeV for all the light pseudoscalar 



mesons of the octet. This is suggested by the sum rule for g which shows that g/ f is 
flavour independent. In eq. ( ^.1| ) we have detailed the contributions A°, A^ and the 



one from Ag/g respectively. We have also taken Aqsb = 1 GeV. It is evident that 
there is a large cancellation among the last term and the other ones. In order to be more 
quantitative we have to fix the value of g, which, on the basis of our result (|3.23|) , depends 



on the value of F. In Ref. [^ the range of values g ~ 0.2 — 0.4 was found; therefore, 
putting f7^ = 0.1, we would obtain 



^B 



1.3 MeV (4.2) 



The application of our results to the charm case is more doubtful, in view of the large 

values of the l/rric correction (a — b)/mc. By scaling the result ( [4.2| ) to the charm case, 

one obtains 

An = —Ab ~ 4.4 MeV (4.3) 

rric 

In conclusion, our estimate of gp*p*T, — gp*p-K allows to include a previously neglected 
term in the loop induced contribution to the hyperfine splitting. Although our estimate 
is affected by an uncertainty in the value of Gs, nevertheless this new term tends to cause 
a substantial cancellation and to reconcile the chiral calculation with the experimental 
data. 
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